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We demonstrate that at a finite temperature, an effective wall thickness of a single walled carbon 
nanotube (SWNT) should he W = Ws + Wd, where Ws is the static thickness defined as the 
extension of the outmost electronic orbit and Wd the dynamic thickness due to thermal vibration 
of atoms. Both molecular simulations and a theoretical analysis show that Wd is proportional to 
\/T. We find that the increase of dynamic thickness with temperature is the main mechanism of 
Raman spectrum shift. The introduction of the dynamic thickness changes some conclusions about 
the Young's modulus and reduces the values of thermal conductivity. 

PACS numbers: 



a 



I 

O 

o 



en 
> 

00 

o 
o 



I 

O 

o 



% 



Nanotubes have attracted increasing attention in the 
last decade due to potential applications in nanoscale 
electronic, mechanical and thermal devices [1-13]. De- 
pending on the geometrical structure, nanotubes can ex- 
hibit fascinating properties, for example, by varying the 
chiral index {n,m), nanotubes can change from semi- 
conductors to metals. Y-K Kwon et al.[lj| studied the 
thermal contraction effect of fuUerenes and carbon nan- 
otube. They find that in nanotubes, the gain in entropy 
translates into a longitudinal contraction, which reaches 
a maximum at 800K. Schelling et al.[l5j studied the 
thermal expansion coefficient of carbon nanotubes with 
empirical bond-order potential. These studies show that 
the structure of carbon nanotube depends strongly on 
temperature. In addition to tube radius and length, wall 
thickness is another important and fundamental struc- 
ture parameter. It can infiuence the Young's modulus 
and thermal conductivity, two typical physical quanti- 
ties characterizing the mechanical and thermal properties 
of nanotubes. These two quantities depend crucially on 
the cross-section of the tubes. Therefore, any ambiguity 
in calculating the thickness will cause an error in these 
quantities, and might results in misleading conclusions. 
Unfortunately, up to now, to our best knowledge, a uni- 
fied and unique way to define the thickness of nanotube 
does not exist even though some discussions have been 
madepBlllllllllillllrnl. In the existing literature, a 
so called static thickness, defined as the extension of the 
outmost electronic orbit, is used as the wall thickness. 
This static thickness is temperature independent. 

In this paper, we will demonstrate, with numerical evi- 
dence and theoretical arguments, that at a finite temper- 
ature the static thickness alone cannot be used as the wall 
thickness of SWNTs. One needs to introduce a dynamic 
thickness to refiect the thermal effect. 



To carry out the measurement of wall thickness, a cam- 
era with resolution on the atomic scale is supposed to be 
put along the axis of SWNT. The wall thickness is deter- 
mined by measuring the width of the wall image. If the 
camera shutter is in femtosecond (fs), a clear picture will 
be obtained, and the wall thickness from this image is the 
instantaneous static thickness; but if the shutter time is 
longer than fs, such as in picosecond (ps) scale, the ef- 
fect of the thermal vibrations of atoms will appear, the 
wall image becomes blurry and the wall appears much 
thicker. We know that in most physical process such as 
heat conduction along the tube[8|,lll|, the time is much 
longer than ps, so the widening of the wall thickness from 
thermal effect can not simply be ignored. We call the con- 
tribution from this part dynamic thickness and denote it 
as Wd ■ The wall thickness should he W = Wg + Wd- 
Thus, in calculating any relevant physical quantity, such 
as Young modulus and the thermal conductivity etc, one 
should use the effective thickness W instead of the static 
thickness Wg- We will demonstrate this with numerical 
simulations that the dynamic thickness is a natural struc- 
ture parameter and affect many physical properties such 
as Raman spectrum shift. Young's modulus and thermal 
conductivity. 

We begin our discussion by studying the radial distri- 
bution of SWNT at different temperatures. The results 
are obtained from molecular dynamics (MD) simulations 
in which the Tersoff empirical bond order potential 18] is 
used. The Hamiltonian of the carbon SWNT is: 
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where Vij — 



fc{rij)[VR{rij) + hjVA{rij)] is the Tersoff 
empirical bond order potential. Vnirij) and Vyi(ry) are 
the repulsive and attractive parts of the potential, re- 
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FIG. 1: (a) The radial distribution of (5,5) SWNT at differ- 
ent temperatures, (b) Thermal expansion of radius for (5,5) 
(solid square) and (10,10) (solid circle) SWNT. The number 
in the figure is the slope of the lines obtained from best fit- 
ting, (c) The dynamic thickness, VKd, versus temperature for 
(5,5) (solid square) and (10,10) (solid circle). The best fit- 
ting gives rise to: Wd = 7 x io-^t"*^^ for (5,5) tube, and 



Wd 



1.6 X 10"^r" 



for (10,10) tube. 



spectively, and fdr) depending on the distance between 
atoms, bij are the so-called bond parameters depend- 
ing on the bounding environment around atoms i and j 
they implicitely contain many-body information. TersofF 
potential has been used to study thermal properties of 
carbon nanotubes successfully. For detailed information, 
see RefHi]. 

In the MD simulations, the carbon atom is treated as a 
mass point of zero size. The carbon atom vibrations con- 
tain two parts: motion parallel and perpendicular to wall. 
The vibrations perpendicular to the tube wall will intro- 
duce a time dependent SWNT radius distribution. In our 
calculations, three armchair- type SWNTs with different 
diameters, (5, 5), (10, 10) and (15, 15) are simulated, each 
contains 500, 1000 and 1500 atoms, respectively. For a 
given temperature one million steps are used for equi- 
libration and over 30 000 simulation steps are used to 
calculate the radial distribution. The time step is 0.8 fs. 



TABLE I: Dynamic thickness (second column) of three types 
of nanotubes at room temperature compared with different 
static thickness used in different models (columns 3-5). 



Tube Wd Ws 

(5,5) 0.115 A 

(10,10) 0.278 A 

(15,15) 0.304 A 



WilWs 

= 0.7A[12] Ws = 1.44A|iri Ws = 3.4A[13j 

16% 8% 3% 

40% 20% 8% 

43% 21% 9% 



In Figure. 1, we show the radial distribution of (5,5) 
SWNT, thermal expansion of the radius and the tem- 
perature dependence of dynamic thickness for (5,5) and 
(10,10) nanotubes, respectively. The radius distribution 
is well represented by a Gaussian distribution. The center 
of the distribution is the average radius of a carbon nan- 
otube at the corresponding temperature, and the width of 
the distribution is defined as the dynamic thickness, Wd, 
of SWNTs. From Figure 1 we can see that the dynamic 
thickness increases with temperature. The reason is that 
dynamic thickness arises from the thermal vibration of 
atoms. As the temperature increases, the amplitude of 
thermal vibration increases. Wd versus temperature T is 
drawn in Figure 1(c) in double logarithmic scale, which 
indicates that Wd = CT". For both (5,5) and (10,10) 
SWNT, a « 0.49, which agrees very well with the follow- 
ing theoretical analysis. 

In a SWNT, each carbon atom has three nearest neigh- 
bor carbon atoms which are bond together by covalent 
bonds. The central carbon atom vibrates perpendicular 
to the plane determined by its three nearest neighbors 
(see Fig. 2). If the central atom deviates from its equilib- 
rium position, the force arising from its nearest neighbors 
will drag it back. Within the first order approximation, 
this vibration can be treated as an harmonic one with 
an effective spring constant k^. If the amplitude of each 
oscillator is Ai?,, then (Ai?^) « iksT/ke, where ks is 
the Boltzmann constant and T temperature, ke is the ef- 
fective spring constant for the thermal vibration and can 
be calculated as we will show below. Thus the average 
width caused by the thermal vibration is 
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The calculated dynamic thickness Wd, for three differ- 
ent nanotubes at room temperature is shown in Table 
I. These values are compared with the static thickness 
of the corresponding nanotube used in previous stud- 
ies [13, ll3 UM (columns 3-5). Dynamic thickness in- 
creases as the radius increases, while the increase de- 
creases for large tubes. This arises from their different 
curvatures. Fig. 2 shows the detailed structures of one 
carbon atom with its three nearest neighbor atoms in 
(5,5) and (10,10) SWNTs. The SWNT can be seen as 
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FIG. 2: Structure sketch for (5,5) and (10,10) SWNTs. 



a wrapped-up graphite sheet. In a graphite layer, four 
carbon atoms share a common plane; while in a SWNT, 
the central carbon atom is not in the same plane with its 
three nearest neighbors. The distance between the atom 
and the plane determined by its three nearest neighbor 
atoms is represented hy h. h changes with tube diam- 
eters and is proportional to the curvature of the tube 
wall, thus, the larger the tube radius the smaller h be- 
comes (see Fig. 2). For (5,5) SWNT, h is about 0.151, 
and for (10,10) SWNT, h is about 0.07A. Each carbon 
atom in a SWNT is considered to be connected with its 
nearest neighbors through the covalent bond K'^^'^. K'^^'^ 
is the effective spring constant for the covalent bond. It 
is determined by the atomic interactions. K'^~'^ is larger 
for strong covalent bond than for weak bond. It can be 
shown approximately that the effective spring constant 
for the thermal vibration is, ke oc K'^~'^h? . Therefore, at 
a given temperature, the SWNTs with larger radius have 
smaller /i, and larger dynamic thickness. 

The \/T comes from the MD simulation results. This 
relation can also be deducted from the simple harmonic 
model we describe above. The perfect consistent between 
the simple harmonic model and the numerical simula- 
tion demonstrate that in the temperature range we study 
(< 500K), the vibration of carbon atom in SWNT is 
harmonic, this is consistent with other theoretical study 
that only when T > SOOKthe anharmonic effect will 
appear[l^. The MD simulations is also performed for 
a zigzag tube and similar effect of dynamic thickness is 
found. So the dynamic thickness is a natural structure 
parameter to all nanotubes. In fact, from the analysis 
above, we can see that the dynamic thickness comes from 
the vibrational of atom, so it exists in all types of nan- 
otubes. 

There is not a unified view on the static thickness. For 
example, 0.617 — 0.77A was used in a continuum mechan- 
ics model |l3;|l6J, 1.44A llTl was used for the diameter of 
carbon atom, and 3.4A llJl, the inter-layer separation of 



graphite, was also used as the wall thickness. These dif- 
ferent thickness are static thickness only, while for effec- 
tive thickness, one should include the dynamic thickness. 
From Table I, one can see that at room temperature, even 
for the smallest dynamic thickness of (5,5) SWNT, the 
dynamic thickness varies between 3% to 16% of the static 
thickness; while for bigger (15,15) SWNT, the dynamic 
thickness can be as high as 10% ~ 40% of the static 
thickness. Moreover, the dependence of thermal proper- 
ties of SWNT on temperature are important problems in 
carbon nanotubes studies|j,|^ |3,l3- Therefore, when a 
temperature-independent static thickness is used to cal- 
culate the thermal conductivity, the information about 
the temperature dependence in thermal conductivity will 
be lost. The high ratio and strong temperature depen- 
dence show that the dynamic thickness can not be ig- 
nored in the wall thickness calculation. Ignorance of the 
dynamic thickness will induce a large error in calcula- 
tions of relevant properties such as Young's modulus and 
thermal conductivity. In the following we shall use the 
increase of the wall thickness due to the thermal effect to 
explain the Raman spectrum shifty a correction of thermal 
conductivity and Young's modulus, respectively. 

Ram,an spectrum, shift. In high-temperature studies 
of Raman-active modes in SWNTs, it is found that the 
Raman peak frequency shifts toward lower frequency as 
temperature is increased |l9[ l2(l |. In an attempt to de- 
scribe this shift quantitatively, a temperature coefficient 
of Raman frequency 
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was introduced 19j. It is observed experimentally that 
a^ = -2.47 X IQ-^K-'^ for a SWNT with diameter of 
1.34nm which is close to a (10,10) SWNt[3|. This fre- 
quency shift has been explained as the diameter ther- 
mal expansion of SWNTs and the softening of the intra- 
tubular bonds [l9j. However, the molecular dynam- 
ics simulation in Ref [l9j gives a much larger value of 
a^{= —5 X 10^^). Here we attribute the Ramman fre- 
quency shift to the thermal expansion of the radius, afj, 
and the increase of the thickness, a^"^ . As we shall see, 
we obtain results which correspond more closely to ex- 
periments. 

After taking into account the dynamic thickness, the 
Raman radial breathing frequency uj goes as l/(i?-f M^/2), 
where R is the radius of tube and W is the wall thick- 
ness of tube. As Wg is independent of temperature, the 
temperature coefficient can be written as: 
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dR 1 dWd 



R + W/2 \dT 2 dT 



a^: + a^^ 



(4) 



Generally speaking, the contribution from thermal ex- 
pansion of the radius, af^, is about one order of magni- 
tude smaller than that one from the increasing thickness 
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FIG. 3: Temperature coefficient of Raman frequency shift 
ttui versus temperature T for nanotube witli radius (5,5) and 
(10,10). 



a^-". For example, for a (10,10) nanotube at T == SOO/sT, 
a^, w -0.1 X IQ-'^'K"^ and a,^'' « -2.84 x IQ-'^R-^, 
thus we have a^^ « —2.94 x lO^^A'^^, -which is close 
to the experiment value, —2.47 x IQ^^K^^ . Our results 
show that the temperature effect of dynamic thickness 
is another source of Raman spectrum shift besides the 
bond softening effect. 

Moreover, as dW/dT ^ 1/VT, the absolute Raman 
frequency shift is larger at lower temperature and become 
smaller at higher temperature as is illustrated in Fig. |3 
where the a^j versus T is shown. We believe that this 
temperature dependence behavior should be observed in 
experiment. 

Thermal conductivity. Thermal conductivity is very 
sensitive to wall thickness, as heat current is defined as 
energy transport through a unit cross-section area in unit 
time. Different thicknesses will give rise to different cross- 
sectional area of the nanotubes, thus different heat cur- 
rent and thermal conductivity. The thermal conductiv- 
ity, corrected after introducing the dynamic thickness, is 



K = Kp 1 



Wd 



Ke (l + /3Vt) , 



(5) 



where Ke is thermal conductivity calculated with the 
static wall thickness Wg and /3 = 2y/2kB/ke/Ws. The 
thermal conductivity. Kg, is calculated from a non- 
equilibrium molecular dynamics method, namely, two 
thermal baths with slightly different temperature are put 
into contact with the two ends of the nanotube, after a 
sufficiently time, a stationary state is reached and a tem- 
perature gradient is set up, the conductivity is then calcu- 
lated by using, k = — JL/AT. Where, J is the heat flux, 
L the tube length, and AT the temperature difference 
of the two heat baths. For more information about the 
calculation, please see Rcf. 9] Thus the real thermal con- 
ductivity is always smaller than the experimental one. In 
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FIG. 4: Thermal conductivity for (5,5) nanotube calculated 
with static thickness {Ws), fi:e (solid circles), and with effec- 
tive thickness {W — Ws + Wd), k (solid triangles), (a) Ther- 
mal conductivities calculated with Ws ~ 0.7A. (b) Thermal 
conductivities calculated with Ws = 1.44A. The tube length 
is fixed at 50 atomic layers. The curves in the figure are drawn 
to guide the eyes. 



the worst case, when Wd/Wg ^ 0.4, the true thermal con- 
ductivity is approximately 70% of the measured value. In 
other words, the thermal conductivity in SWNTs calcu- 
lated by using static thickness, is somehow exaggerated. 
In Fig. ^we show the thermal conductivity k and k^ for 
(5,5) nanotubes with different static thicknesses. Fig^ 
(a) corresponds to thermal conductivities calculated with 
Ws = 0.7 A while Fig|3|(b) with Ws = 1.441. These cal- 
culations clearly demonstrate that the correction of the 
dynamic thickness is very significant, in particular in high 
temperature regime. 

Young's modulus. Another important property of nan- 
otubes is the Young's modulus which determines its me- 
chanical property. In Yakobson's calculation with (7,7) 
SWNT O, W^s = 0.7l was used, the Young's modu- 
lus calculated was 5.5 Tpa; while Lu[lJ| obtained 0.97 
TPa with Ws = 3.4A. The difference can be removed if 
they use the same wall thickness, thus these two Young's 
modulus are in fact consistent with each other. One shall 
note that if we only discuss the Young's modulus of the 
same tube, the wall thickness is not very important, but 
if we compare the Young's modulus of different tubes, 
the effect of dynamic wall thickness will be very impor- 



tant. For example, in Ref. [13, by using Wg = 3.4A for 
both (5,5) and (10,10) SWNTs, with empirical force con- 
stant model calculation, the Young's modulus for (5,5) 
and (10,10) SWNT are 0.971TPa and 0.972TPa, respec- 
tively. From this result, it is concluded that the Young's 
modulus is insensitive to the radius of the SWNT. How- 
ever, if we notice that at the same temperature, the wall 
thickness is not uniform for (5,5) and (10,10) SWNT, 
then the conclusion is different. By using the dynamic 
thickness of both tubes at room temperature (see table 
I), and Ws = 3AA as the static thickness, one obtains the 
Young's modulus, 0.939TPa and 0.899 TPa for (5,5) and 
(10,10), respectively The smaller SWNT has a slightly 
larger Young's modulus than the larger SWNT. If the 
smaller Ws is used, then the dynamic thickness correc- 
tion would be even more obvious. 

We would like to point out that, in experiments, it is 
impossible to measure directly the thickness of SWNT. 
The thickness is indirectly inferred from the measure- 
ment of other properties, such as the Ramman frequency. 
The wall thickness is an "effective thickness" , namely the 
dynamic thickness is already included. More recently, a 
hardness transition was observed by changing pressure on 
nanotubeJ21j. In this study, it is obtained, from room- 
temperature MD simulations, that an effective wall thick- 
ness is about 0.66A, which should already include finite 
temperature effect. In the current paper, we split the 
effective wall thickness into static thickness and dynamic 
thickness so as to give the wall thickness a clear physics 
picture. 

In summary, we have introduced a dynamic thickness, 
Wd, to SWNTs. Unlike static thickness, which is inde- 
pendent of temperature, the dynamic thickness is tem- 
perature dependent and increases with temperature, T, 
as vT. At room temperature, Wd is comparable to the 
carbon atom diameter, and thus cannot be ignored in 
calculating the thermal conductivity or the Young's mod- 
ulus for different tubes and/or at different temperature. 
The introduction of dynamic thickness can alter previous 
conclusions wherein only static thickness is used. More- 
over, we have found that the increase of dynamic thick- 
ness with temperature is the main mechanism of the Ra- 
man spectrum shift. Our numerical calculations agree 
very well with the experiment data. The dynamic thick- 
ness has more significant effects when the static thick- 



ness is small, such as in SWNTs. If the static thickness 
is large, the effects comes from dynamic thickness will 
decrease. And for certain measured properties, the ef- 
fects of dynamic thickness will be weakened because of 
the time averaged in measurements. We believe that dy- 
namic thickness is a natural structure parameter with 
real physical effect for nano scale systems and can be 
measured experimentally such as in the Ramman fre- 
quency shift. 
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